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Abstract 

The Tutte polynomial of a graph is a 2- variable polynomial which is quite important 
in both combinatorics and statistical physics. It contains various numerical invari- 
ants and polynomial invariants ,such as the number of spanning trees, the number 
of spanning forests , the number of acyclic orientations , the reliability polyno- 
mial, chromatic polynomial and flow polynomial . In this paper,we study and gain 
recursive formulas for the Tutte polynomial of pseudofractal scale-free web(PSW) 
which implies logarithmic complexity algorithm is obtained to calculate the Tutte 
polynomial of PSW although it is NP-hard for general graph. We also obtain the 
rigorous solution for the the number of spanning trees of PSW by solving the recur- 
rence relations derived from Tutte polynomial , which give an alternative approach 
for explicitly determining the number of spanning trees of PSW. Further more,we 
analysis the all-terminal reliability of PSW and compare the results with that of 
Sierpinski gasket which has the same number of nodes and edges with PSW. In 
contrast with the well-known conclusion that scale-free networks are more robust 
against removal of nodes than homogeneous networks (e.g., exponential networks 
and regular networks). Our results show that Sierpinski gasket (which is a regular 
network) are more robust against random edge failures than PSW (which is a scale- 
free network) Whether it is true for any regular networks and scale-free networks 
,is still a unresolved problem. 
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1 Introduction 



The Tutte polynomial is a two-variable polynomial which can be associated 
with a graph, a matrix, or, more generally, with a matroid . This polynomial 
was introduced by W.T. Tutte p~H3] and has many interesting applications 
in several areas of sciences such as combinatorics, probability, statistical me- 
chanics and computer science. It is quite interesting since several combina- 
torial, enumerative and algebraic properties of the graph can be investigated 
by considering special evaluations of it [3J. For instance, one gets information 
about the number of spanning trees [5j[6] , spanning connected subgraphs [7J, 
spanning forests [8] and acyclic orientations [9] of the graph by evaluating 
tutte polynomial at particular points (x, y) . Moreover, the Tutte polynomial 
contains several other polynomial invariants, such as flow polynomial [JTI] , 
reliability polynomial [TT] and chromatic polynomial [2JIT2JH2]. It has also 
many interesting connections with statistical mechanical models such as the 
Potts model [131416] and the percolation [17] . Despite its ubiquitythere are no 
widely-available effective computational tools to compute the Tutte polynomi- 
als of a general graph of reasonable size. It is shown that, many of the relevant 
coefficients do not even have good randomised approximation schemes and var- 
ious decision problems based on the coefficients are NP-hard [T81IT9] .Although 
it is hard to to compute the Tutte polynomials, a lot of efforts have been de- 
voted to the study the Tutte polynomials of different graph such as polygon 
chain graphs [IB], Sierpinski gaskets J2U] and strips of lattices [2"IH2"3] .However, 
polygon chain , lattices and Sierpinski gaskets cannot well mimic the real-life 
networks, which have been recently found to synchronously exhibit two strik- 
ing properties: scale- free behavior [25] and small- world effects [26] . 

Pseudofractal scale-free web(PSW) we studied is a deterministically growing 
network introduced by S.N. Dorogovtsev [27] which is used to model scale-free 
network with small-world effect. Lots of job was devoted to study its properties 
,such as degree distribution , degree correlation , clustering coefficient [271 128] 
, diameter [28]. average path length [29], the number of spanning trees |30j and 
mean first-passage time for random walk [31j.As for its tutte polynomial and 
reliability polynomial ,to the best of our knowledge, related research was rarely 
reported . 

In this paper,we study and gain recursive formulas for the Tutte polynomial 
of PSW .The analytic method is based on the its recursive construction and 
self-similar structure .Recursive formulas for various invariants of Tutte poly- 
nomial can also obtained based on their connections with the Tutte polynomial 
. We also obtain the rigorous solution for the the number of spanning trees by 
solving the recurrence relations derived from the Tutte polynomial,which coin- 
cides with the result obtained in [30] .Thus we give an alternative approach for 
explicitly determining the number of spanning trees of PSW .Further more,we 
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analysis the all-terminal reliability of PSW and compare the results with that 
of Sierpinski gasket which has the same number of nodes and edges with 
PSW. In contrast with the well-known conclusion that scale-free networks are 
more robust against random node failures than homogeneous networks (e.g., 
exponential networks and regular networks) .Our results show that Sierpin- 
ski gasket (which is a regular network) are more robust against random edge 
failures than PSW (which is a scale-free network) . 



2 Preliminaries 

2. 1 The Tutte polynomial 

Let G = (V(G), E(G)) denotes a graph with vertex set V(G) and edge set 
E(G); we will often write V and E, when there is no risk of confusion, and 
so G = (V, E). A subgraph H = (V(H),E(H)) of a graph G = (V(G),E(G)) 
is said spanning if the condition V(H) = V(G) is satisfied. In particular, a 
spanning tree of G is a spanning subgraph of G which is a tree. Let if be a 
spanning subgraph of G and k(H) be the number of connected components of 
if. then the rank r(H) and the nullity n(H) of H are defined as 

r{H) = \V{H)\ - k{H) = \V(H)\ - k{H) 

n{H) = \E(H)\ - r{H) = \E(H)\ - \V(H)\ + k(H) 

Definition 1 Let G = (V, E) be a graph. The Tutte polynomial T(G; x, y) of 
G is defined as ^ 

T(G;x,y) = £ (x - i)'(G)-r(/0 (y _ iy(H) (1) 

HCG 

where the sum runs over all the spanning subgraphs H of G . 

The Tutte polynomial can be evaluated at particular points (x, y) to give nu- 
merical invariants, such as the number of spanning trees, the number of forests 
and the number of connected spanning subgraphs. The following theorem [1] 
depicts the relations between the Tutte polynomial and its numerical invari- 
ants. 

Theorem 1 Let G = (V, E) be a connected graph , Then: 

(1) T(G; 1, 1) is the number of spanning trees; 

(2) T(G; 1,2) is the number of spanning connected subgraphs of G; 
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A 

Fig. 1. Growth process for PSW from n=0 to n=2 
(3) T(G; 2, 1) is the number of spanning forests of G; 

The Tutte polynomial also has a variety of single-variable polynomial in- 
variants associated with the graph,such as the all-terminal reliability poly- 
nomial,flow polynomial and the chromatic polynomial. In this paper we only 
study the all-terminal reliability polynomial. Let G = (V, E) be a graph,each 
edge of G has a known probability p of being operational; otherwise it is failed. 
Operations of different edges are statistically independent, while the nodes of 
G never fail. The all-terminal reliability polynomial R(G,p) of G is defined as 
is the probability that there is a path of operational edges between any pair of 
vertices of G.The connection between the Tutte polynomial and all-terminal 
reliability polynomial is given by the following theorem |3]. 

Theorem 2 Let G = (V, E) be a graph, Then 

R(G,p) = p^-\l - p )\E(G)HV(G)\ + l T 1; _1_^J (2) 

2.2 Structure of PSW 

The scale-free network we studied is a deterministically growing network which 
can be constructed iteratively [27] . We denote the pseudofractal scale-free 
web(PSW) after n iterations by G(n) with n > 0. Then it is constructed as 
follows: For n = 0, Go is a triangle. For n > 1, G n is obtained from G n -\\ 
every existing edge in introduces a new node connected to both ends of 
the edge. The construction process of the first three generation is shown in 
Fig. [TJ The PSW exhibits some typical properties of real networks. Its degree 
distribution P(k) obeys a power law P(k) ~ ^i+in3/in2 ^ e aV erage path 
length scales logarithmically with network order [29|. The network also has 
an equivalent construction method [30,33j,as can be seen in Fig. [2]: to obtain 
G(n + 1), one can make three copies of G(n) and join them at the three most 
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Fig. 2. Second construction method of PSW: G(n+1) is composed of three copies of 
G(n) denoted as Gi,G 2 ,G3 ,the three hub nodes of which are represented by A,B,C 
in the corresponding triangle(The left side of the figure ).In the merging process, hub 
node A of G\ and hub node B of G3 ,hub node A of G3 and hub node B of G2 , hub 
node A of G2 and hub node B of G\ are identified as a hub node A,B,C of G(n + 1) 
respectively (The right side of the figure ). 

connected nodes denoted by A,B,C, which are called the hub nodes of G(n+1) 
in this paper. According to the second construction algorithm, one can see that 
at each step n, the total number of edges in the network increases by a factor 
of 3. Thus, the total number of edges for G(n) is E n = 3" +1 .We can also find 
that the total number of nodes for G{n) is V n = 3 " 2 +3 . 



3 Tutte polynomial of PSW 



Let us simply denote by T n (x,y) the Tutte polynomial of PSW G(n) .In 
this section we study and gain recursive formulas for T n (x,y).The analytic 
method is based on the the relation between spanning subgraphs of G(n + 1) 
and spanning subgraphs of G(n). According to the second construction algo- 
rithm,we find that there exists a bijection between spanning subgraphs of 
G{n + 1) and spanning subgraphs of G\, G 2 , G 3 inside G{n + 1), while the sub- 
graphs G±, G2, G3 are isomorphic to G(n). Indeed, if H is a spanning subgraph 
of G(n + 1) ,we can uniquely determines three spanning subgraphs Hi, H2 
and if 3 of G±, G2 and G3 ; viceversa, given three spanning subgraphs Hi, H 2 
and H 3 of Gi, G2 and G3, respectively, then their union provides a spanning 
subgraph H of the whole G(n+1). Therefore, according to the Definition [TJthe 
Tutte polynomial of PSW G(n + 1) can be rewritten as 

T n+ i(x,y)= £ (x-l)«^'(!/-ir (if) , (3) 

HiCGi,i=l,2,3 



where Hi is the spanning subgraph of Gi ,and H is the union of Hi, if 2 and 
H3.I11 order to obtain the recursive formulas for T n (x, y), we want to know the 
relations between r(G(n+ 1)) and r(G(n)) ,r(H) and r(Hj), n(H) and n{Hj), 
for % = 1,2,3. It is easy to know that r(G(n+l)) = 3r(G(n))-l and \V(H)\ = 
\V(H X )\ + \V{H 2 )\ + \V(H 3 )\ - 3, \E{H)\ = \E{H X )\ + \E{H 2 )\ + \E(H 3 )\ ,for 
every spanning subgraph H of G(n + 1). Furthermore, two possibilities can 
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occur. 

1) In the spanning subgraph Hi of Gi, the hub nodes of G(n + 1) belong to the 
same connected component, for any i, 1 < i < 3, then 

k(H) = A;(^i)+A;(ff2)+fc(fl3)-2 and r(A) = r(ifi)+r(fT 2 )+r(if 3 )-l 
Thus 

ra(if) = |£(#)| - r(H) = n(ifi) + ra(if 2 ) + n(# 3 ) + 1 

r(G(n + 1)) - r(#) = £(r(G(n)) - r(^)) 

i=i 

Hence, one gets: 

[X - iy(G(n+l))-r(H)^ y _ ^n(H) 
3 

= (v - 1) n( x - i) r ( G ( n »- r ( H % - i) n(Hi) (4) 

i=l 

2) For certain z, 1 < i < 3 , in the spanning subgraph H iy the hub nodes of 
G(n + 1) do not belong to the same connected component, we have 

k(H) = k(H 1 ) + k(H 2 ) + k(H 3 )-3 and r(A) = r(H x ) + r(H 2 ) + r(H 3 ) 

Moreover 

n(H)=n(H 1 )+n(H 2 )+n(H 3 ) 

Hence 

r ( G („ + 1)) - r(tf) = E(r(G(n)) - r(#)) - 1 

i=i 

Thus 

( x _ l)»'(G(n+l))-r(fO(j / _ 

= f\(x - l)r(G(n))-rm( y _ l)n(^) (5) 

fa - 1) f=l 

Let D n denotes the set of spanning subgraphs of G(n),we define the following 
partition on D n : 

• Di n denotes the set of spanning subgraphs of G(n), where the three hub 
nodes belong to the same connected component; 

• .D^n denotes the set of spanning subgraphs of G(n), where the hub nodes 
A and B belong to the same connected component, and C belongs to a 
different component. Similarly, (or D® n ) denotes the set of spanning 
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subgraphs of G(n), where A (or B )does not belong to the connected com- 
ponent containing the other two hub nodes; 
• D 3n denotes the set of spanning subgraphs of G(n), where the three hub 
nodes belong to three different connected components. 

Thus, for any n > 0, we have 

D n = D ljn U D* n U Dl n U D c 2 n U D 3 , n (6) 

For n > 0,let us define polynomial of x, y 

ri,„(x,y)= £ (x- l)Kr»)-r(H) (y _ 

Similarlywe can also define polynomials :T 2 A n (x, y), T 2 B n (x, y), T 2n {x, y), T 3jn (x, y), 
while the sum is conducted on D^, D^, D 2n , D 3iU respectively. For symme- 
try, we have 

and we can simply use T 2>n (x,y) to denote one of the three polynomials. Ac- 
cording to Definition [1] , we have: 

T n (x, y) = T ljn (x, y) + 3T 2>n (x, y) + T 3>n (x, y) (7) 

Furthermore, we obtain the following theorem based on the relation between 
spanning subgraphs of G(n + 1) and spanning subgraphs of G(n). 

Theorem 3 For n>0, the Tutte polynomial T n (x,y) of G(n) is given by 

T n (x, y) = T hn (x, y) + 3(x - l)P n (x, y) + (x - l) 2 Q n (x, y) (8) 
where Ti tn (x,y), P n (x,y),Q n (x,y) satisfy the following recurrence relation: 

T 1>n+1 (x, y) = (y - 1)T 3 „ + 3(y - l)(x - l)7? >n P„ 

+3(x - l) 2 (y - l)T hn N 2 n + (x- l) 3 (y - 1)P 3 
+3(x - l)Tl n Q n + 6Tf >n P n + 12(x - l)T 1)n P n 2 
+6(x - 1) 2 P 3 + 3(x - l) 3 P 2 Q n + 6(x - l) 2 T hn P n Q n (9) 



P„+i(x,y) =4T lin P 2 + 4(x - l)T hn P n Q n + 4(x - 1)P 3 

+4(x - 1) 2 P 2 Q„ + (x - l) 2 Ti,„Q 2 + (x - l) 3 P n Q 2 (10) 
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Q n+1 (x, y) = 8P 3 + 12(x - l)P 2 Q n + 6(x - l) 2 P n Q 2 n + (x - 1) 3 Q 3 (11) 

with initial conditions 

Ti, (x, y) = y + 2 Po(x,y) = Qo(x,y) = I. 

and Ti tn ,P n ,Q n is shorthand ofTi tTl (x, y), P n (x,y) ,Q n (x,y) respectively. 

Proof: We find that x — I divides T 2tn (x, y) and (x — l) 2 divides T 3jn (x, y) for 
any n > which will be proved later. As a consequence, we can write 

T 2 , n (x,y) = {x- l)P n (x,y) 
T 3>n (x,y) = (x-l) 2 Q n (x,y) (12) 

where P n (x,y) and Q n (x,y) are polynomials of x,y. Thus we obtain Eq.® 
from Eq.(j7j). 

Now we will proof Ti tn (x,y), r 2)n (x, y),T 2jn (x, y) satisfy the following recur- 
rence relations which lead to the results of the theorem. 

7\, n+1 (x, y) = (y- l)2* n + 3(y - l)T 2 n T 2 ,„ + 3(y - l)T 1)ft 7f jn 

+(y - l)T 3 n + -^(3T 2 n T 3 , n + 6T 2 n T 2)ri 
x — 1 

+12T ljn T 2 2 n + 6T| n + 3T 3 , n T 2 n + QT^ n T 2)n T^ n ) (13) 
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2,„ +1 (x, y) = —-^T^T^ + 4T lin T 2in T 3 , n + 4T 3 , 



3,nJ 



(14) 



T 3 , n+1 (x, y) = — ^ (8T 3 n + 12T 2 n T 3 , n + 6T 2 , n T 3 2 n + T 3 3 n ) (15) 

with initial conditions 

Ti, (x,y) = y + 2 T 2j0 (x,y) = x - 1 T 3i0 (x, y) = (x - l) 2 . 

and Ti in ,T 2jn ,T 3jn is shorthand of T 1>n (x,y), T 2 , n (x, y),T 3>n (x, y) respectively. 

The strategy of the proof is to study all the possible configurations of spanning 
subgraphs i/j in the three copies of G(n) inside G(n + 1), for i = 1,2,3, 
and analyze which kind of contribution they give to T l n+1 (x, y), T 2 n+1 (x,y) 
and T 3>n+ i(x,y). 



8 




Fig. 3. The possible configurations of spanning subgraphs H{(i = 1,2,3) for 
Ti jn+ i(x, y).The two hub nodes of Gi are connected by a solid line if they are in in 
the same connected component, and connected by a dotted line if they are not in 
the same connected component. 




Fig. 4. The possible configurations of spanning subgraphs Hi{i = 1,2,3) for 
T 2n i i(x, y).The two hub nodes of Gi are connected by a solid line if they are in in 
the same connected component, and connected by a dotted line if they are not in 
the same connected component. 

For Ti n+ i(x, y),we find it has 10 possible configurations of spanning subgraphs 
Hi( i = 1,2,3) which is shown in Fig. EJln the first configuration, Hi G D l n 
for any i = 1,2,3. This contributes to T^ n+ i(x,y) by a term (y — l)Tf n 
according to Eq.(@|, since in the spanning subgraph Hi of Gi, the hub nodes 
of G n+ i belong to the same connected component,for any i, 1 < i < 3. In the 
second configuration, if, G -D l n holds for two i G {1, 2, 3} and Hi G D 2n holds 
for the last z(for example,^ G Di n ,H 2 G Di n ,H 3 G D^ n ).This contributes 
to Ti jn+1 (x, y) by a term 3(y — l)T 1 2 n T 2i „ according to Eq.(j3J. Computing the 
contributions to T ljn+1 (x, y) of the 10 possible configurations and adding them 
together ,we obtain Eq. ffTB"]) . 

For r 2i „ + i(x, y), we study T 2 n+1 (x, y) only by symmetry. We find it has 6 possi- 
ble configurations which is shown in Fig.HlIn the first configuration, H 3 G Di >n 
Hi G Dj n U D^ n ,H 2 G D^ n U D 2n . This contributes to T^ n+l (x, y) by a term 
^jT 1>n T2 n according to Eq.(j5]), since in the spanning subgraph H 1 or H 2 , 
the hub nodes of G n+ \ do not belong to the same connected component. In the 
second configuration,!^ G D l n Hi G D^ n U D 2n ,H 2 G _D 3in .This contributes 
to T^ n+1 (x,y) by a term -^Ti in T 2jn T 3:n according to Eq.flHJ). Computing the 
contributions to T 2n+l (x, y) of all the 6 possible configurations and adding 
them together ,we obtain Eq.flHJ). 

For T^^ n+ i(x, y), we find it has 4 possible configurations which is shown in 
Fig. [5jln the first configuration, Hi G D 2 n U D 2 n ,for any i = 1,2, 3. This 
contributes to T 3 n+1 (x, y) by a term -^T 2n according to Eq.(j5]), since in the 
spanning subgraph Hi H 2 and H 3 , the hub nodes of G n+ i do not belong to 
the same connected component. In the second configuration,^ G D 2n U 
holds for two i G {1,2,3} and Hi G D 3n holds for the last z(for example, Hi G 
Di n U D$ H 2 G D$ n U Df tn ,H 3 G A,,„).This contributes to r 3>n+1 (x,y) by 
a term ^-T 2 2 n T 3 n . Computing the contributions to T 3n+1 (x, y) of all the 4 
possible configurations and adding them together ,we obtain Eq. (jT5l) . 
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Fig. 5. The possible configurations of spanning subgraphs H{(i = 1,2,3) for 
T3 jn+ i(x, y).The two hub nodes of G{ are connected by a solid line if they are in in 
the same connected component, and connected by a dotted line if they are not in 
the same connected component. 

For the initial conditions, It is easy to verify according to the definition. 

Now, we come back to proof the recurrence relations that Eqs. (f9|) . (TTQl) . (fTTT) 
show. First,we find x — 1 divides T 2jn (a;,y) and (x — l) 2 divides T 3in (x, y),for 
n > 0. The results can be proved by mathematical induction based on the 
recurrence relations that Eqs. (fT4^ . ffT5]) show. 

Inserting Eq.([12]) into Eqs. (USD,®, (USD for T 2jTl (x,y) and T 3>n (x,y),we obtain 
Eqs.fl9]). ffT0^ f lTTj) .The initial conditions for P n (x,y) and Q n ( x ,y) is easy to 
verified according to the initial conditions for T 2>n (x, y) and T 3 n (x, y). 

Remark:Although it is NP-hard to calculate the Tutte polynomials for gen- 
eral graph,the recurrence relations we obtain shows that we can calculate the 
Tutte polynomials for PSW with time complexity 0{n) = O (log {V n )). Thus 
,we have obtain logarithmic complexity algorithm to calculate the Tutte poly- 
nomial of PSW . 

We can also obtain the recursive formulas for various invariants of Tutte poly- 
nomial based on their connections with the Tutte polynomial, such as the num- 
ber of spanning trees,the number of connected spanning subgraphs, the number 
of spanning forests,the number of acyclic orientations ,the reliability polyno- 
mial and the chromatic polynomial .In this paper ,we only study the number 
of spanning trees and the reliability polynomial. 



4 Exact result for spanning trees of PSW 

Let us denote by Nsrin) the number of spanning trees of PSW G{n) .Accord- 
ing to Theorem [1] and Theorem [3j it is easy to know 



Further more , for any n > 0,we have the following recurrence relation 



N ST {n) = T n (l,l) =Ti >n (l,l) 



N ST (n+l) = 6N ST (n) 2 P 7 



n 



and 



P, 



n+l — 



4iV ST (n)P n 2 



where P n is abbreviation of P n (l, 1) ,and the initial conditions is 



iV sr (0)=Ti,o(l,l) = 3 



1. 
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Thus 

N ST (n) = QN ST (n-l) 2 P n „ 1 

= 6^ 1+2 h 1 [Ns T (n-2)f x2+ ^Pi 2 _ + 1 2) 

^6 a H bk [N ST {n- fc)] c *P*i 
where a>k,bk,Ck,dk,k > 1 satisfy the following recurrence relations 

ttfc = dk-l + Cfc~l 

bk = bk-i + dfc-i 
c fc = 2cfc_i + 4_i 
4 — c fc-i + 2c4-i 

with initial conditions 

di = 1 61 = Ci = 2 and di = 1 

Thus 

c k + d k = 3(c fc _x + 4_i) = 3 fc_1 (ci + di) = 3 fc 
c k - d k = c k -i - d k -i —d — di — 1 

Then 

3 fc + 1 



Cfc = 



dk = 



2 

3 fe - 1 



Note 

afc - a-k-i = Ck-i and 6 fc - b k _ x = d k -i 

Thus 



fc-i fc-i k — 1 3 k — 3 
a fc - ai = XX a *+i - «i) = c * = ~~ ^ 1 7 — 

i=l i=l 

fc-i fc-i fc — 1 3 k — 3 
bk-h = J2( b i+i -bi) = J2 d i = 7— + — ] — 

i=i i=i z 4 
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Hence 



k + l , 3 fc -3 

«fc = + (25) 

fc-1 3 fc -3 

fefc = + — ^ (26) 



Substituting Eqs.(|23|>,(|23D, ([2SD,(I2SDinto Eq.dH]),we obtain the following re- 
sult. 

Theorem 4 For any n > 0,the number of spanning trees of PSW G(n) is 
given by 

, , 3 n+l_ 2n _ 3 3 n+l +2n + 1 

N ST (n) = 2 3 3 3 (27) 



Note:The result coincides with the result obtained in [30], and we give an al- 
ternative approach for explicitly determining the number of spanning trees 
for PSW. It is smaller than the the number of spanning trees for Sierpinski 
gasket [32] . 



5 Reliability analysis of PSW 

5.1 All-terminal reliability of PSW 



In this section, we look upon PSW G(n) as a probabilistic graph. Each edge 
of G(n) has a known probability p of being operational; otherwise it is failed. 
Operations of different edges are statistically independent, while the nodes of 
G(n) never fail. The all-terminal reliability polynomial R(G(n),p) of G(n) is 
defined as is the probability that there is a path of operational edges between 
any pair of vertices of G(n). In general case ,the calculation of all-terminal 
reliability polynomial is NP-hard [31]. But for PSW, we obtain recursive for- 
mulas for all-terminal reliability polynomial which implies that logarithmic 
complexity algorithm is obtained. Further more,we get a approximate solu- 
tion of R(G(n),p) based on the recursive formulas it satisfy, which shows that 
all terminal reliability decreases approximately as a exponential function of 
network order. 



let us simply denote by T l n ,P„,the expression T l ri (l, jz^),P n f 1 5 yz^J respec- 
tively for n > 0, we obtain the following recurrence relations from Theorem 
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3. 

r -( 1 'rb) = ri -( 1 «rh;) (28) 

U = ^ + 6T 1 2 n P n (29) 

^n+l f 1 ' l3^J = 4T l,« P n (30) 

with initial conditions 

Jl^]^ Pod, 1 ^ • 



1 — p J 1 — p \ 1 — p J 

For n > 0,let us simply denote by i?(n) the reliability polynomial R(G(n),p) 
of PSW G(n) ,and define 

B(n) = p (v ^ 2) (l - pf E -- v - +2) P n (31) 
We obtain the following recurrence relation from Theorem [2] 



R(n + 1) =p v '"+i- 1 (l - p) E "+ 1 - y «+ 1 + 1 T 1)n+1 (l ' 



1 — p 

p 3V„-4 (l _ p) 3 En - 3 y„ + 4 ( _^_ T 3 n + 6T 2 nPn) 



,3(V„-1) 



+6p 2(V„-l) (1 _ p) 2(^-V tt+ l)^ . p (V„-2) (1 _ p) {E n -V n+ 2)p n 

R(n) 3 + 6R(n) 2 B(n) (32) 



and 



B(n + 1) =jo (y " +1 - 2) (l - p )(E n+1 -v n+ i+2) Pn+i 

= 4p (3^-5) (1 _ p) (3E„-3K 1+ 5) Tijn p2 

= 4i?(n)5(n) 2 (33) 

where V n ,E n denote the total number of nodes and edges of G{n) respectively 
and the initial conditions is 

R(0) =p 2 (3-2p) B(0) =p(l -p) 2 

Note that R(n) ^> B{n) while n — > oo,ones get 



13 



R(n + 1) + 2B(n + 1) = R(n) 3 + 6R(n) 2 B(n) + AR(n)B(n) 2 

= (R(n) + 2B(n)f - A(R(n) + 2B(n))B(n) 2 
^(R(n) + 2B(n)f (34) 



Thus 

R(n) « (i?(n - 1) + 25(n - l)) 3 

^( J R(1) + 25(1)) 3 "- 1 = [pp-p)] 3 "" 1 « [p(2-p)]f y " (35) 

which shows that all terminal reliability decreases approximately as a exponen- 
tial function of network order V^.The reason we don't use R(n— l) 3 as a approx- 
imation of R(n) is that it has lager relative error than (R(n — 1) + 2B(n — l)) 3 . 
In fact, we find that B(n) is the probability that G(n) is split into two different 
connected components such that one of them contains the hub nodes A,B ,the 
other one contains the hub node C .Thus R(n) + 2B{n) < l,for any n > 0. 

5.2 Comparison of all-terminal reliability between Sierpinski gasket and PSW 

The Sierpinski gasket is a fractal which can be constructed iteratively [55] . 
The starting point is a triangle. Divide its three sides in two segments of equal 
length. Connect the midpoints to get four inner triangles and paint the three 
external ones. Apply the same process to the inner triangles but the middle 
one. Sierpinski gasket is the limiting set for this construction. 

If we look upon the Sierpinski gasket as a network,it is a deterministically 
growing network which has the same starting point with PSW .But the method 
of iteration is different from PSW |36j. The construction process of the first 
three generation is shown in Figj6jWe denote the Sierpinski gasket after n 
iterations by SG(n) with n > O.It has the same number of nodes and edges 
with PSW for any n > 0,but the structure is quite different. For Sierpinski 
gasket, except the 3 outmost nodes which have degree 2, all other vertices of 
SG(n) have degree 4. In the large n limit, SG(n) is 4-regular.But PSW is a 
scale free netwokr whose degree distribution obeys power law .Thus, Sierpinski 
gasket and PSW are typical examples of regular network and scale free network 
which have the same number of nodes and edges . 



A A 

SG{0) SG(1) SG(2) 
Fig. 6. Growth process for Sierpinski gasket SG(n) from n=0 to n=2 . 
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Alfredo Donno [20J found that for each n > 0, the all-terminal reliability 
polynomial of Sierpinski gasket SG(n) is given by 

R(SG(n),p) =p v "- 1 (l-p) E "- v " +1 T hn (l, 1 



1 — p 

with T n (l, j^) = T 1>n (l, j^) and 

Ti,n+i ( 1, -^-1 = r^-^n + 6I? >n JV„ (36) 
\ 1-pJ l-p 

N n+ i ( 1, ^- J = T^-^n + 3t„Af n + 7T hn N 2 (37) 
\ 1 — p j 1 — p ' 

M n+1 f 1, — !— ) = -^T lin iV n 2 + 12T ljn iV n M n + 14iV n 3 (38) 
with initial conditions 

For n > 0,let us simply denote by R s (n) the reliability polynomial R(SG(n),p) 
,and define 

B s (n) = p {v ^ 2 \l - pY E "~ v " + ^N n (39) 
T s (n) = p^-^il - p)^ En ~ Vn+ ^M n (40) 

We obtain the following recurrence relation . 

R s (n + 1) =p v ^-\l - p) E ^- v ^ +1 T 1>n+l (l, — !— ) 

1 — p 

= p3 y„-4 (1 _ p) 3^- 3 y„ + 4 ( ^ T 3 + 6T 2 Nn) 

1 — p ' 

= R s {n) 3 + 6R s {n) 2 B s {n) (41) 



and 



B s {n+l)=p (yn+1 - 2 \l _ p )(^ + i-v„ +1 +2) iV) 



n+l 



:p (3Vn-B) (l _ p) (3E^3V n+5){ _P_ T 2 nNn + + 77^) 

: R 8 {nfB s {n) + i?» 2 T» + 7R s (n)B s (n) 2 (42) 
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T s (n+ l)=p {Vn+1 - 3) (l - p) iEn+1 - Vn+1+3) M n+1 

= p {3V n -6)^ _ p )(3B n -3V n +6)( JP-T \ >n Nl + 12T hn N n M n + UN 3 ) 

1 — p 



= 3R s (n)B s (n) 2 + 12 R s (n)B s (n)T s (n) + UB s (n) 3 (43) 



with initial conditions 



R s (0) = p 2 {3-2p) B a (0) =p{l-pf T s (0) = (l-p) 3 



Now, we compare all-terminal reliability between Sierpinski gasket and PSW 
while p G (0, 1) . For n=0, 



Thus R s (l) = -R(l), But for n > 0,we can obtain from Eqs. (I33|) and (j42p that 
B s (n) > B(n). Thus ,for any n > l,we find from Eqs. fl32l )and fj41H 



We have calculated R s (n) and R(n) for different p G (0, 1) and n by itera- 
tion. we find that R s {n) and R(n) converge to quickly with n — > oo and 
Eq. (l44p holds for n > l.The results for n = 6 are shown in Fig. [TlOur 
results show that Sierpinski gasket is more robust than PSW against ran- 
dom edge failures. Thus,we obtain an example which shows that regular net- 
works^. g., Sierpinski gaskets) are more robust than scale-free networks(e.g. 
,PSW) against random edge failures. Whether it is true for any regular net- 
works and scale-free networks ,is still a unresolved problem. However, recent 
work [3"TH3"U] have shown that inhomogeneous networks, such as scale-free 
networks, are more robust than homogeneous networks (e.g., exponential net- 
works and regular networks) with respect to random deletion of nodes. Thus, 
combining with our above result, we can reach the following conclusion that 
networks (e.g., scale-free networks) which are more robust again random node 
failures do not mean more robust again random breakdown of edges than those 
(e.g., regular lattices) which are more vulnerable to random node failures . 



6 Conclusion 

In this paper,we study and gain recursive formulas for the Tutte polynomial 
of PSW which implies that recursive formulas for various invariants of Tutte 
polynomial can also obtained based on their connections with the Tutte poly- 
nomial . We also obtain the rigorous solution for the the number of spanning 



R a (0) = R(0) 



B s (0) = B(0) 



R a (n) > R(n) 



(44) 
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Fig. 7. The all-terminal reliability of Sierpinski gasket SG(6) and PSW £7(6) ob- 
tained by direct calculation from Eqs. ([52"j) and (T4T]) . 

trees of PSW by solving the recurrence relations derived from Tutte polyno- 
mial , which give an alternative approach for explicitly determining the number 
of spanning trees of PSW. Further more, we analysis the all-terminal reliability 
of PSW based on the the recurrence relations derived from Tutte polyno- 
mial and compare the result with that of Sierpinski gasket. In contrast with 
the well-known conclusion that scale-free networks are more robust than ho- 
mogeneous networks (e.g., exponential networks and regular networks) with 
respect to random deletion of nodes. Our results show that there is an example 
that regular networks(e.g., Sierpinski gaskets) are more robust than scale-free 
networks(e.g.,PSW) against random edge failures Whether it is true for any 
regular networks and scale-free networks ,is still a unresolved problem. 
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